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EFFECT OF STABILIZING FORCES ON TURBULENCE* 
By L. Prandtl 



I have frequentljr noticed on pleasant summer evenings that 
a temperature layer, such as is caused hy the evening radiation 
and the ensuing cooling off of the air strata near to the 
earth's surface, manifests itself in the sudden cessation of 
the turhulence of the wind and flows on without any. sign of 
forming clouds. There is "but one interpretation for this, 
namely, the warm air glides in laminas over the cold layers of 
air "beneath it. The appearance of expressed discontinuity 
layers in the free atmosphere, where a warm air mass flows over 
a cold mass without perceptihle intermingling, "belongs in the 
same group, so also the flow of fresh water over sea water as 
is particularly observed in the polar ocean. The opposite 
phenomenon is that of increased turbulence due to convective 
motions ensuing from weight differences when the earth and the 
"bordering lowest layers of air or water are warmed up. In the 
first case the "exchange" as well as the friction of the flow- 
ing air or water mass is lowered, in the second case it is 
raised. Since these interrogatory forms were in close relation 
to our &Bttingen test program on fully' developed turbulence,** 
I acquiesced to the request of Wilhelm Schmidt of Vienna to 
undertake a study of this subject. In these experiments an 
air stream is to be blown jDetween a water-cooled and a vapor- 
heated plate. According to what preceded it is anticipated 
that the exchange will be lowered appreciably or even stopped 
altogether under appropriate conditions if the stratum is 
staple (warm above, cold below), while in the contrary case 
the exchange will be enhanced.. Because these experiments had 
to be made at relatively low air speeds so as to let the lift- 
ing forces become effective, it first required a study of the 
test methods, and so we have nO;t.yet proceeded beyond the pre- 
liminary tests on the most suitable hot-wire connections and 
thermocouples for recording speed and temperature. The wish 
to predict the phenomena which we want to observe here prompt- 
ed me, however, to search into the theory of these phenomena. 



*"Einfluss stabilisierender Kraft e= auf die Turbulenz , " -.f rom a. 
reprint of Vortrage aus dem Gebiete der Aerodynamilc und 
verwandter Gebiete, Aachen, 1929, ..issued by A. Gilles, L. Hopf 
and Th. v. Karman. 

**N.A.C,A. Technical Memorandum lio. 435. Sept.. 1926, ."Turbulent 
Flow." 
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In addition I was hopeful that this aspect vrould he productive 
of enhanced insight into the mechanism of turbulence. 

The momentary status of this theory is as follows: 
Assume a horizontal fluid flow of vertical laminas so that each 
superposed stratum is lighter than the one "below it. Then, 
when the motion is turbulent it produces energy, e.g., the 
heavier layer is raised and the lighter layer is lowered against 
the lift of the heavier. Let the path in the vertical, -traversed 
"by a fluid portion prior to intermingling anew with its new sur- 
rounding be designated by I, ' This is the so-called "mixing path 
or distance." If -dp/dy expresses the decrease in density in 
the vertical (y- direct ion) , then the ultimate lift difference 

of a particle shifted upward becomes- = "^^^y ' -.Since this lift 
difference Increases from zero to this value while covering dis- 
tance I, we have an energy outpti.t of *^^S^^fy (positive when, 

as assumed her'e, dp/dy is negative!). Downward motions yield 
the same ez]pression. In order to embrace the volumes partici- 
pating on the exchange, we draw a horizontal surface, and at- 
tribute to one fraction thereof an upward motion with speed 

Vj^, and to another fraction Ps a downward motion with mean 
speed V2. Accordinci;ly the -total amount on surface F is 
^(piVi + psv^) (volume per second) , thus making available as per- 
formance the product of the "apparent shearing stresses" called 
forth by the exchange and the rate of slippage. When we bring 
this performance into relation with a body of base P and 
height l (Pig. 1 ) , 

pTt 

- , : — _ — > 

Fig-1 



and when T» is the apparent shear stress and u the mean rate 
of the flow, this performance becomes = ?T'*l ~. The selection 

of a layer of height I is. .justified, for the previously con- 
ceded fluid quantity is precisely transm.itted in such a layer. 

Conformal to 0. Eeyrxolds the apparent shear stress is = pu'^v'^, 
the line denoting a formed average. Accordingly the mean of 
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u' = and of v« is = 61 v, + pgVg. Since the sign for v' 

dy 

fluctuates regularly with that of u' we write 



The existence or nonexistence of energy for the further 
maintenance of the flow in question now depends on which of the 
two energy figures predominates. If the energy output exceeds 
the shear stress, the difference is available for upholding the 
turbulence, hut if the lift performance is greater than the 
weight difference, the turbulence must perforce die away. Thus 
something is left to maintain the turbulence when 

FT' > r(p^v^ + P3V3) t6^"dy> , 

after inserting the value of T' and aVor eviat ing 

dy. 2 dy ■ 

One notable feature is the elimination of the chosen sur- 
face F, the mixing distance I and the mean value of verti- 
cal speed from our inequation (l), so that this interpre- 
tation does not reveal anything about these two important 
quantities. In addition, it is- advi sable to designate the 
quotient at the right and left sides of (l) separately. Making 
this nondimensional quantity =0: 

e = — — ^y-, . (2) 

yields t^irbulence for 6< 1 and strong stability; the turbu- 
lence dies out when. 6 > 1. G =0. denotes that the fluid is 
homogeneous. Sitice'the derivation of the preceding formulas 
is primarily a rough est.im.ate, the experiments undoubtedly call 
for some kind of correction factor? But. from, the construction 
of the formulas- it may be* inferred ' that ; for instance, what the 
relevant quantities are and what qualitative course may.be ex- 
pected. An inherently homogeneous warm air stream flows over 
an identically cold one; for instance, du/ dy may temporarily 
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te very high, in the transition layer .iDptvreen the two, so as to 
comply with inequation (l), that" is, prevailing turbulence. 
During the turhulent -exchange the mixing distance "becomes, grad- 
ually thicker, and du/ dy conformably smaller. This also ap- 
plies to the transition of the temperature, Uow the subsequent 
limiting attitude; is plausible ;;s?here inequation (l) becomes pre- 
cisely the equation and then protracts the attained attitude be- 
cause of the discont inxied turbulence (the pure viscosity effects 
and the thermal conductivity are to be small enough to become 
insignificant).-. .N.ow, with h as the height of the transition 
layer in. ?;?hdch .ll;Ae^^ distribution prevails, we write 

du • -"^s ."^1 1 do 

*^ and" instead of ^ ± ^iL with T = absolute tern- 



• ■ • :. h- p dy 

p erasure, we. put i 4^ =: -L. ^S. 

T dy Tj^ h 

This yields 



or 



2(u2 - Ui)^T 



m 



g(T2 - 



(3) 



This 'equation (s) gives the minimum value of the height of such 
a transition layer at which the flow is without turbulence. The 
comparison with empirical figures should show the accuracy of 
the numerical factor of this formula' or any eventual correction, 
although a precursory collation attested to the conformity of 
(3) with experience. 

The shear stress T» of the apparent friction nay be ex-* 
pressed by 

T' = pi^(|3i)%(e), (4) 

dy 



regarding which we, for the time being, merely say that f(e) 
must become =^ 0 when 0^1, and f(0) must become = 1 
when 9=0 (ordinary turbulence in homogeneous fluids). On 
the premises of an energy- consideration (intensity of turbulent 
motion proper t lo.nal to the difference, of the above explained 
energies), It is -.concluded that f(e) A^r.~;^*g, subject how- 

*The difference of the potential t emperatures'^shouldT^'s t ^ 
speaking, be used here' for this temperature difference so as to 
embrace the terapernture changes due to changes in heights as 
the pressure varies. 
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ever to further proof. Ahove all it is anticipated that nixing 
distance I will te affected ty 9, a thought which might "bo 
extended to the case of initially unstable stratification, "or 
in other words, that the air "below is warmer than above. Such 
conditions may arise for lon^^cr 'periods when the air is warmed 
up from the ground* In this case the mixing distances will 
presumably be greater due to the brisk thermal convection than 
in the case of ordinary turbulence. 



Closely bound up with these conditions of laminar flow are 
the curved flows, of a homogeneous fluid, in our particular case 
the plane flow in concentric circles. Let the mean speed in 
tangential direction u' be a function of the radius. When the ■ 
speed decreases from the outside toward the inside the fluid 
particles from the inner zone can no longer penetrate the outer 
flow, for their centrifugal force is less than that- of the outer 
particles, and conversely can no longer hold the particles which 
came from without, but return whence they came. Conversely, 
instability prevails when the inward speed rises sufficiently. 
To each particle motivated frictionless in the remaining medium 
the surface theorem expressed as u = tangential component of the 
speed u r = constant is applicable,' because the pressure field 
is precisely symmetrical in rotation. As soon as a particle is 
radially shifted for an amount dr., its speed obviously changes 
by - '^dr. If the mean speed at the new position differs from 



r 

old 



the old by ^dr, a radial slippage for an amount I yields in 
dr ■ 

first approximation a speed difference 



u 



dr 



r 



(5) 



This difference becom.es zero for the potential flow with circu- 
lation, u = and consequently plays the same role with 

r 

respect t o • tu.rbulent motion as the motion with constant speed in 
the rectilinear flow. A particle which differs in speed- 
from its surrounding shows a difference in centrifugal 
unit volume with respect to that in the vicinity 



u» 
force 



per 



p(u + uO 



- u^ 



is, in first , approximation, 



that 

. r 

Thus we find in connection wi.th formula .(5) that potential . 
motion u r = constant is an exact parallel of the case of 
neutral . equil ibrium. An outward increase of p'r.oduct • ur-. de- 
notes. a stable arra^ngdment, an outward : decrease , an unstable 
process. "In the lattdr case we encounter, the vortex motions ' 



s 
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of which G. I, Taylor had made a particular study,* 

How we revert to the two energy output s . As far as the 
energy against /the centrifugal force is concerned this is for 
a slippage of I = l/2 times the defined difference in cen- 
trifugal force per unit of voliime, or 

plUU» 2U /dU ^ Us 

The amount of the additive fluid is computed as before. When 
defining the adduced output it must "be r eniemher ed that the 
speed of the relative slippage of two intersecting surfaces 
(speed of deformation) distant hy the length unit is here 

d^ ^ (equal zero for the rotation as. rigid hody u = a;r) . 

Then the output "becomes P tH (~ - , where again 

dr r 

T' = pu'^v"^' - pUg + |) (j3^v^ -f P2^3-^' 

and the formula, expressing the possibility of turbulence 
finally .becomes 

(d:U ^ U) (du ^ uy ^ (du H. u ^ (g) 
dr r dr r dr r r 

If ^.+ 5 i3 positive, it yields |5 > 2^, if M + ^ 
dr r - " _ ^ dr r 

= negative, we have ~ < 2-. But the latter inequation, 

d.r r _ 

positive .u assumed, is always complied with for ~ < - ^ 

dr r 

according to as s.umpt Ion . 

The conditions are more conveniently perceived when 
graphing the right and left sides of ineq\iation (5), first 

dividing by' (— ) • and. using ~ : — "as ,absc issa , as trans- 
r ■ • cir r 

scribed on Figiire 2, where the 'zones of anticipated turbv.lence, 
i.e., where (6) is fulfilled, are shown as shaded portions* An 
examination of temporary typical . floyys whose speed distribution 
conforms to power law 's: .=■ ar^ yields:. 



♦Taylor, G.. I..: "Stability, of a Viscous . Liquid Contained between 
two Ro-tating . Cylinders , " Phil. Trans. Roy. Soc, Vol. 223 A, 
1923, p. 28S. 
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instability,. Taylor vortices for n < - 1 

indifference, potential flow " n = - 1 

stability, laminar flow - 1 < n < + 2 

stability, rigid rotation " n + 1 



stability , 



n > + ^ 




Fig. 2 



naturally, it nuist be proved bj^ experiment whether the 
limit for the turbulence in' the^ outward increasing speeds lies 
actually at n = 2 or at sone other figure. However, precur- 
sory studies revealed that- the- limi.t- value is not far from 2. 

In order to obtain a: fornula for the shear stress we 
again have recourse to the ratio of- the right side of inequa- 
tion (6) to the left side as parameter, so that when making 
it = we can write 

Then fCG^) = 0 for ^ 1 and = ^"^"^7 ^^r smaller 

values.* The latter reference Is yet without tangible in- 
formation so long as we do not know- how the mixing length I 
is dependent on 6,, It; follows that the aim should be to so 
define function ffg^^.) that I would be contingent on the 

geometrical" configuration "of the cylinder but not on the speed 
distribution. The possible ambit- of this hinges on the studies 
under preparation to this end. 

* It becomes readily apparent that G^^ = 0 is conformal to the 

rectilinear motion (r = ^).. so that our formula here reduces 
again to the usual form. . • . 
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As regards these experiments they pertain, in a few words, 
to water hetween two concentric cylinders, each of which can he 
driven at" a different rate of rotation. The torque transmitted 
through the water is to he measured and the speed and. pressure 
distrihution defined. 



Commentary hy J. Mi Burgers, Delft 



The effect of a density layer on the stability of flow was 
studied by L. F. Richardson.* His deductions relate to the 
prohlem of atmospheric flow with thermodynamic relations and he 
adduces the decay of the turbulence when 

0 dh Cp Bh ah • 

Here 6 = absolute temperature, g = acceleration of gravity, 
Cp = specific heat of air at constant pressure, h = height, 

v^j. and Vy = the horizontal components of the mean speed. 

With 0=9+ 9^ becomes the potential temperature 

and the first term of the above inequation assumes the form 
g i& p ' * ■ . ■ . 

Q "^^^ A series of observations made, by Benson revealed a 

satisfactory affirmation of the formula. It is manifest that 
this criterion agrees with Prandtl's reference'. Hichardson^s 
studies may be applied direct to a simple case where the 
thermal expansion is without any significance. 

I may be allowed to add a few remarks relative to the 
formula for the exchange coefficient, particularly when the 
mean flo^ has only one speed component u in direction x 
which is a function of y, all other quantities in question 
to be functions of y. According to definition we then, ob- 
tain the value of the exchange coefficient in the neighborhood 
of a point P by laying a surface S through P parallel to 
the x-y plane and defining at moment t for each point of 
this surface the momentary value of v' (i.e., the speed com- 
ponent in direction y) and of distance I* traversed by the 
volume element passing through the relevant point since the 
"last compensation process" in positive direction y. These 
distances P may take a different turn according to whether 



*L. F. Richardson: Proc. Roy. Soc. London, Vol. 97, 1920, p. 
554. Phil. Mag. Vol. 49, 1925, p. 81. 
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we stress the compensation of the mean speed, or temperature, 
or of the concentration of a detached substance, etc. Each 
special case then yields the pertinent exchange coefficient as 
the mean value 



over surface S. 

Now the expression "last compensation" is not easily de- 
fined, hut we may, if we first study the compensation of the 
•speed, arrive at a satisfactory formula as indicated hy 
G. I. Taylor,* when we substitute for I' the distances I 
which the relevant volume elements have traversed since a 
stated time interval T(<t) in the positive direction y, 
Hereshy interval t T must "be chosen large enough to ensure 
at least one. speed compensation in each volume element. 
Whereas v* and I' (i.e., the speed in direction y at 
moment t and the distance traversed since the last compensa- 
tion) always have the same sign, there must, if the conception 
"last compensation" is correctly ixaterpreted, exist no correla- 
tion "between V and I - I*. Consequently 

vW" - = v"T - v"r~-"T"y = v"T. 

The last quantity is unob j ec t ionably defined and may, according 
to Taylor, be used in the formulas of the speed exchange. 

Taylor, moreover,, related quantity v' I with two other 
quantities. It is patent that 



(because forming. the mean value always pertains to the attitude 
at. moment t on the fixed surface S, thus permitting the 
change in consecutive order of differentiation and mean values). 

On the other hand 

•' ^^'t""" t 

V»T = (V^-^ / V^T-dT) =7 V»^v«T.dT, 

T " T- 



♦Taylor, G.I. : Proc. London Math. Soc, Vol. 2, 1922, p. 196. 
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when v*^ is the speed in direction y , at interval T, of the 

same volume element to.T7hich v*^ is relevant; (the formation 

of the mean value is always with respect to surface S in time 
interval t ) . 



With 



we have 



t t«-T 

V»l = V»2 / R^^T*^'^ ^ / ^1^^' 



is a so-called correlation coefficient, v/here the fact that 
R^ is a linear fxmction of has "been utilized. It naturallj^ 

is now assumed that / Rfcd^ has a finite value; it is really 

0 ^ 

this limiting value which defines the value of v'l. 
Consequently 

^ dt 0 ^ 

Quantity is related with a similar one occurring in 

the theory of BrownVs motion.: On the other hand, mean value 

"^t^T "belongs to the type of "correlation moments" introduced 
"by Keller and Priedmann in their .general theory on turhulence.* 

Heretofore we were concerned with the exchange of speed 
or momentum. How we turn to the exchange of an arbitrary 
quantity and express the relevant exchange coefficient hy 

Now we must distinguish two cases: the exchange of the 
considered quantities is slower than that of the speed or that 
of the relevant quant it ies - is faster than that of the speed. 
In the former, l\ is general ly measur ed from an earlier time 

interval than I*, thus precluding any correlation "between 
V » and l\ - I ' • Then 



*Eeller, L. and A. Friedmann: Proc. 1st. Intern. Congress 
Applied Mechanics, Delft, 1924, p. 395. 
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In the latter case, l^^ is measured from a later time interval 

than I', so that I' - usually has the same prefix as v», 

and we have: 

€^ . = - v"rT~-T»iy < v"T^ or < G. 

The derivation of the Richardson formula is hound up with 
the exchange coefficient for density gradients. The volume 
elements passing in interval t through S have the density 

and are, in consequence, pulled downward with a force 

""''i dh " 

The force acting against the gravity in unit time per unit 
of volume is 

- gv^l 



i dh 



So when we presume the density gradients to compensate at least 
not faster than those of the speed, we can substitute € for 

With respect. .to the last point it should he borne in mind 
that the coopensation of differences in density by pure diffu.sion 
is usually slower than that of- speed differences by internal 
friction. Because K/\ic^ is greater than 1 for many gases, 

the compensation of temperature differences may proceed more 
rapidly than that of the speed differences by internal friction 
alone; but, since the latter compensation is .aided by the pres- 
sure effects no danger, is likely to occur. 



Translation J. Vanier, 
National Advisory Committee 
for Aeronautics. 



